Equivalence between focussed paraxial beams and the quantum harmonic oscillator 
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The paraxial approximation to the scalar Helmholtz equation is shown to be equivalent to the 
Schrodinger equation for a quantum harmonic oscillator. This equivalence maps the Gouy-phase of 
classical wave optics onto the time coordinate of the quantum harmonic oscillator and also helps us 
understand the qualitative behavior of the field and intensity distributions of focused optical beams 
in terms of the amplitude and probability distributions of quantum harmonic oscillators and vice 
versa. 
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I. INTRODUCTION 

It is well known that in the limit of short wavelengths, 
wave optics can be reduced to ray-optics: Maxwell's 
equations can be reduced to the eikonal equations of ge- 
ometric optics. |3| Ray optics therefore is an emergent 
theory with limited applicability. The ray optical descrip- 
tion of optical image formation is simple and satisfactory 
for many purposes. Indeed, if thin lens formulae apply 
and complications such as chromatic aberration can be 
neglected, the description is so simple that it is standard 
material in undergraduate physics education. But op- 
tics requires a wave description when the characteristic 
length scales of a beam of light become comparable to its 
wavelength. Even when the thin lens formula appears to 
be applicable, a wave optical description becomes neces- 
sary where a focused light beam's transverse extension 
shrinks to a point, for example, in all focal areas of a 
conventional imaging apparatus. Wave optics limits the 
resolution of point-images through aperture dependent 
point-spread functions. 0, 13, S IE 

We investigate the fundamental building blocks of gen- 
eral images: focused, monochromatic beams. Specifi- 
cally, we consider focused beams in the paraxial approx- 
imation and assume that we can treat the polarization 
degrees of freedom separately so that we can apply the 
paraxial approximation to the scalar Helmholtz equation. 
Since a general image is described by a suitable mixture 
of coherent beams (of different frequencies, polarization, 
etc.) our considerations are of some generality. Beams 
that are not too strongly focused can typically be de- 
composed into a complete orthonormal set of modes that 
have a simple analytical description. The most widely 
used are the Hermite-Gaussian transverse electromag- 
netic modes (TEM- modes) , which we will introduce in 
Sec. m These modes are essentially suitably rescaled 
wave functions of a harmonic oscillator multiplied by 
phase factors arising from the specific geometry of the 
beam, see Eq. Q.lillll This connection between co- 
herent optical beams and harmonic oscillators is known, 
but it seems to be little appreciated that the two sys- 
tems are isomorphic (but see Ref. |3)- The amplitude of 
a focused paraxial beam corresponds to the amplitude of 
a quantum harmonic oscillator, the former's intensity to 



the probability distribution of the latter, and the Gouy- 
phase of optics 8] assumes the meaning of time for the 
quantum harmonic oscillator. 

We derive this equivalence by a transformation of the 
paraxial wave equation of optics into the Schrodinger 
equation for the quantum harmonic oscillator in Sec. lIIII 
We then give a few illustrative examples for this equiv- 
alence in Sec. IIVI which will allow us to discuss some 
counterintuitive aspects of optical beam behavior. 



II. HERMITE-GAUSSIAN BEAMS: 
TEM-MODES 

General monochromatic TEM-beam modes in the 
paraxial approximation are solutions of the associated 
scalar Helmholtz equation for the electro-magnetic vec- 
tor potential A propagating in vacuum at speed c (with 
the dispersion relation, uj = cfc)(Q] 
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where r = (x, y, z) and t describe position and time 
coordinates. A monochromatic beam with wave num- 
ber k travelling in the positive z-direction and lin- 
early a;-polarized is described by the vector potential 
A = (Ax,Ay,Az) whose only non-zero component is A^ 
with0| 



A,(r,t;fc) =^(r) e'C^'-'^^l 



(2) 



If we insert this ansatz into the Helmholtz equation 
and factor out the common plane wave phase factor 
gi(fez-wt)^ we obtain a differential equation for the field 
envelope ip{r) alone: 
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Equation Ij3b|l represents the paraxial approximation; it 
is based on the observation that, for beams that are not 
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too strongly focused (several wavelengths across) , the z- 
dependence of the field derivative is primarily due to the 
plane wave factor 2iA:, which allows us to discard the 
dependence on dtl^/dz. When we perform the second 
derivative, the term d'^tp/dz'^ is consequently assumed to 
be negligible. 

The paraxial approximation significantly simplifies our 
task because it has the familiar transverse electromag- 
netic TEMmn-modes Vmn as solutions. 0, 0, 0| They 
contain products of Gaussians and Hermite polynomi- 
als in the transverse beam coordinates x and y, that is, 
the familiar harmonic oscillator wave functions ipm {£,) = 
(Oexp( " 



if™(Oexp(-eV2)/x/2^ 
various phase factors. 
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The relation that links the minimal beam radius wq with 
the Rayleigh range b via the light 's wa vele ngth A and 
wave number k is given by wo = ^/2b/k = ^/Xb/ir. The 
beam radius (where the intensity has dropped to 1/e^ of 
the central intensity) at the distance z from the beam 
waist obeys w{z) = ^ywQ{^ + z^/b'^), and for large z im- 
plies the expected amplitude decay of a free wave oc l/\z\ 
with a far-field opening angle arctan(A/(7ru'o)). 

The corresponding wave front curvature is described 
by the radius R{z) — (z^ + b'^)/z, and the longitudinal 
Gouy-phase shift (5|,|^] follows 



(z) = arctan(z/&) , 



(5) 



which varies most strongly at the beam's focus. Note 
that the associated Gouy-phase factor e~*('"+"-+i)'^(^) de- 
pends on the order n and m of the mode functions. Dif- 
ferent modes therefore show relative dephasing or mode- 
dispersion, particularly near the focus. 

The vector potential of Eq. Q describes a beam 
travelling in the positive z-direction (h = kz) and yields 
an electric field that is polarized in the a;-direction with 
a small contribution in the z-direction due to the tilt of 
wave fronts off the beam axis. 6] According to Maxwell's 
equations in the paraxial approximation, that is, neglect- 
ing transverse derivatives, we find for the electric field 
vector,!^ (x,y, z are the unit-vectors and Re stands for 
the real part) 

E(r, t; k) = Re{ [xwV(r) + iic^^]e'^''''^'^ } . (6) 

The associated instantaneous electrical intensity distri- 
bution is la 



/(r,i) = |E(r,t)2 



(7) 



III. TRANSFORMATION TO HARMONIC 
OSCILLATOR SCHRODINGER EQUATION 

We have derived the paraxial approximation of the 
wave equation in Eq. (O and found that the solutions, 
except for some phase factors and a coordinate rescaling 
in the arguments of the solutions, look like those of the 
quantum harmonic oscillators. We shall therefore now 
"undo" these various factors and transform the paraxial 
wave equation into Schrodinger's equation to establish 
the isomorphism between both systems. 
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FIG. 1: Plot of new coordinate lines ^ 
xz-plane. 



constant in the 



Based on the shape of the solutions Q , we apply the 
following coordinate transformations [cc, ?;,z] i— > [^,?7, t], 
which establishes new coordinate lines ^ and i] that follow 
the focused beams' hyperbolic flow lines (see Fig. 
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V2' '^'-V2 
Because of Eq. Q , we let 

z — 6tan(r), 



w{z) 



(8) 



(9) 



where r plays the role of time. If we use 
these new coordinates and the ansatz ^'{x, y, z) = 

^(Cj /w{z) in the paraxial wave equa- 

tion 0, we find after a lengthy but straightforward cal- 
culation 



vl/(e,r7;r) =0, (10) 



which is the Schrodinger equation for the harmonic oscil- 
lator. The units are fi = m = f2 = 1, where m is the mass 
and SI — 27r/T is the resonance frequency of the harmonic 
oscillator. We have established the equivalence between 
focused beams and harmonic quantum oscillators. 

Note that the z-coordinate z = (— oo, oo) in the op- 
tical case was transformed to time r = (— 7r/2,7r/2) in 
the harmonic oscillator case, whereas the role of time in 
the optical case has no analog in the harmonic oscillator 
case because we had to factor out the plane wave phase 
factor e'('^^~'^*). In other words, the isomorphism we es- 
tablished links the beam's envelope amplitude ip{x,y,z) 
with a two-dimensional harmonic oscillator wave function 
*(C,ry;r). 




FIG. 2: Top: Contour plot of the time- averaged intensity 
l{x,0,z) of a field with mode structure TEM40, namely, 
ip4,{V^x /w{z))ipo{\/2y /w{z)) , around the beam focus (0,0,0); 
Rayleigh length b = 20A. In this and all other figures of op- 
tical beams the intensity is in arbitrary units and the coordi- 
nates are in units of A. Bottom: Contour plot of the evolution 
of the probability distribution of a one-dimensional quantum 
harmonic oscillator with wave function 'I'(^) = 'fiiiO- 



IV. EXAMPLES: SINGLE-MODE AND 
TWO-MODE PROFILES 

For simplicity, we will now consider one-dimensional 
beam patterns. We choose the y-mode to be the zero- 
mode, and thus concentrate on the behavior in the trans- 
verse a;-direction only, that is, we investigate the beams 
in the a;z-plane. For easy comparison, we depict focused 
beams with a beam parameter of & = 20A. Our first ex- 
ample is that of a sing le TEM4o-beam, plotted in Fig. |21 
we next consider a more interesting case, a two-mode 
superposition of TEM- modes 30 and 40 (see Fig. 

For the two-mode superposition in Fig. |2| we can 
clearly see the interference between modes 953 and ip^ 
at work and how it leads to the redistribution of the 
intensity from one beam edge to the other, thus in- 
verting the transverse beam profile about the beam 
axis. This intensity redistribution is ultimately due 
to the mode-dispersive effect of the Gouy-phase factor 
^i{m+n+i)<j) -(-j-^g optical case and the time phase factor 
gi(m-i-i/2+n+i/2)nt f^j. ^j^g two-dimcnsioual harmonic os- 
cillator. Note that a beam that is focused only in one 
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FIG. 3: Average intensity I{x, 0, z) of a field with mode struc- 
ture [ip4{^/2x/wo) — itp3{V2x/wo)](po{V2y /wo) at the beam 
focus z = 0; b = 20A. Bottom: Contour plot of the evolu- 
tion of the non-normalized probability distribution of a one- 
dimensional quantum harmonic oscillator with wave function 



transversal direction (using a cylindrical lens rather than 
a spherical one) corresponds to a one-dimensional har- 
monic oscillator, because in this case the Gouy-phase 
factor is e'('"+^/^)'^ and corresponds to a one-dimensional 
harmonic oscillator phase factor g*(™+i/2)f2t^ Clearly, 
both systems evolve through half an oscillation period 
(t = -T/4, . . . , T/4. (Remember, the Gouy-phase © 
varies from — 7r/2 to 7r/2.) 

We know from ray optics that images become inverted 
at the focus. Here we see how this image inversion arises 
through the action of the Gouy-phase. In the language 
of the harmonic oscillator, it corresponds to half an os- 
cillation [— 7r/2,7r/2] or swinging to the "other" side. 

To gain further insight, we plot the beam's instanta- 
neous intensity distribution l{x,0,z;0) in Fig. 0| It can 
clearly be seen that the relative phase of 90° between 
modes 4 and 3 tilts the wave fronts, and the action of the 
Gouy-phase "transports" the beam's envelope across the 
focus. 

Similarly to the case discussed in Fig. |3| but at first 
seemingly less intuitively, we can also have a scenario 
where the harmonic oscillator starts out with an average 
position at the beam center, then swings up to one side 
and falls back to the center. In wave optical terms this 
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FIG. 4: Instantaneous intensity J(2;,0, 2;0) of a field with 
mode structure [(p4{\^x / wo) — i(p3{\/2x /wo)](po{\/2y /wq) at 
the beam focus z = 0. For clarity, a smaller beam parameter 
of b = 5A was chosen. 

scenario corresponds to a field with a symmetric far-field 
distribution that is asymmetric at the focus and becomes 
symmetric in the far-field on the other side, this is illus- 
trated in Fig. |5l 

It also is possible to create a beam that, in oscillator 
language, starts out in the center and swings to either 
side to fall back to the center forming an annulus with 
a dark center at the beam's focus (see Fig. 1^1 . It also 
can be viewed as a beam that is overly focused in the far 
field. Then, according to Heisenberg's uncertainty prin- 
ciple A{x) A{px) > fi/2, such a spatially strongly local- 
ized harmonic oscillator carries a large momentum un- 
certainty. The bottom of Fig. displays the harmonic 
oscillator ground state together with the superposition 
of modes and 2. Clearly the superposition is spatially 
more confined than the ground state, and its Fourier 
transform reveals a large spread, namely, a central peak 
and two momentum peaks that are displaced from zero. 
In other words, this superposition describes a harmonic 
oscillator that "explodes" sideways and then swings back 
to the center. The seemingly counterintuitive formation 
of a dark center at the focus becomes easy to understand. 
Its two-dimensional analog is the optical bottle beam ob- 
served very recently. ^9] The latter is formed from a su- 
perposition of Laguerre-Gaussian modes and 2, and 
thus carries basically the same structure in a cylindrical 
geometry as our example of Fig. |^in the a;z-plane. 

A possible technical implementation of the patterns 
discussed here is sketched in F ig. [71 and has recently been 
demonstrated experimentally, [l^ For further details on 
experimental implementations using computer-generated 
holograms, consult Refs. H m [ij. The field profiles 




FIG. 5: Average intensity I{x, 0, z) of a field with mode struc- 
ture [(p4{\/2x/'Wo) + 'P3{V2x/wo)]'Po{V2y/wo) at the beam fo- 
cus z = 0; b = 20A. Bottom: Contour plot of the evolution of 
the non-normalized probability distribution of a quantum har- 
monic oscillator with wave function \['(^;0) = ip4{(,) + V3(C)- 



sketched here should be useful for optical manipulation 
techniques as well. ^2 



V. CONCLUSION 

Focused beams cannot be described by ray optics in 
their focal region. Therefore an intuitive understanding 
of the wave optical behavior of a coherent, focused beam 
with transverse intensity modulation is desirable. It is 
shown that the paraxial wave equation is equivalent to 
the Schrodinger equation of a two-dimensional harmonic 
oscillator. This equivalence is used to show how the pas- 
sage of a beam through its focus can be understood in 
terms of the behavior of a harmonic oscillator evolving 
through half a period. 

In particular, the inversion of a focused beam's far- 
field intensity pattern about the beam axis, upon passage 
through the focus, is described in terms of a harmonic os- 
cillator swinging to the "other" side, and the formation 
of a bottle beam is explained by the behavior of a con- 
stricted beam pattern that corresponds to a harmonic 
oscillator fanning outward and recoUapsing. 
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FIG. 7: Sketch of possible implementation using a laser, two 
lenses Li and L2, and a suitable hologram H. The hologram 
modulation of the laser beam is mapped into the focal region 
F of lens 1/2- Formally, the transverse field profile in its focal 
plane is the Fourier-transform of the hologram-pattern. 
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FIG. 6: Average intensity I{x, 0, z) of a field with mode struc- 
ture [ipo{V2x/wo) + (p2{V2x/'wo)](po{V2y/wo) at the beam fo- 
cus 2 = 0; 6 = 20A. Bottom: Contour plot of the evolution of 
the non-normalized probability distribution of a quantum har- 
monic oscillator with wave function 'I'(^;0) = 'y5o(C) + '^^iO- 
For comparison this distribution is overlayed by the non- 
normalized distribution for the ground state wave function 
ifoiO alone. It shows that the superposition state ^ is very 
tightly confined in the ^-direction, see text. 
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